Black hole mass dynamics and renormalization group evolution 
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We examine the real-time dynamics of a system of one or more black holes interacting with 
long wavelength gravitational fields. We find that the (classical) renormalizability of the effec- 
tive field theory that describes this system necessitates the introduction of a time dependent mass 
counterterm, and consequently the mass parameter must be promoted to a dynamical degree of 
freedom. To track the time evolution of this dynamical mass, we compute the expectation value of 
the energy-momentum tensor within the in-in formalism, and fix the time dependence by imposing 
energy-momentum conservation. Mass renormalization induces logarithmic ultraviolet divergences 
at quadratic order in the gravitational coupling, leading to a new time-dependent renormalization 
group (RG) equation for the mass parameter. We solve this RG equation and use the result to pre- 
dict heretofore unknown high order logarithms in the energy distribution of gravitational radiation 
emitted from the system. 



INTRODUCTION 

Understanding the dynamics of a system of gravita- 
tionally bound classical or quantum black holes is a prob- 
lem of both theoretical and observational interest. To 
have any hope of finding analytic solutions to this non- 
linear problem necessitates some sort of approximation 
scheme. An obvious avenue for finding such solutions 
lies in the limit where the typical frequency uj of probes 
of this system (e.g. outgoing radiation) is much smaller 
than the inverse of the size of the system and its gravi- 
tational radius r s . 

As we discussed in previous work P, [i| , in this limit 
the dynamics is summarized by an effective worldline La- 
grangian coupled to bulk gravity, describing the evolution 
of a suitable center of mass coordinate x fJ -(X). The world- 
line theory is of the form 



qIt(X)M(X) 



dr(X)l^{X)E l] (X) 



(1) 



In this equation A is an affine parameter, r is proper 
time measured using the bulk metric g^ v and is a 
projection of the bulk Weyl tensor onto its electric-type 
parity components, E fIl/ (X) = C tialJ j3(x(X))x a x l3 . The 
objects M(A), hj(X) describe the first two dynamical £ = 
0, 2 SO(3) multipole moments of the system. Ellipses 
denote higher order moments which do not play a role 
in our present discussion. See 0-Q for a more detailed 
discussion of this Lagrangian. 

The Lagrangian in Eq. ([1]) is suitable for describing ei- 
ther the quantum or classical dynamics of an ensemble of 
black holes in the long wavelength limit. In the quantum 
theory, the moments are operators acting on the Hilbert 
space of black hole microstates. We showed in ref. 
that even without complete knowledge of the microscopic 
theory, it is possible to obtain some information on the 
various correlators of these worldline operators through a 
matching procedure to bulk graviton scattering and ab- 



sorption processes, in a way that is similar to the CFT 
description of bulk black holes in AdS/CFT. 

Alternatively, Eq. ([1]) may be viewed classically, in 
which case the moments are fixed by first treating the 
constituents as point particles and then integrating out 
all of the modes whose wavelength is of order the or- 
bital separation to generate a theory of a single compos- 
ite object. The dynamics of the moments are dictated 
by solving the equations of motion for the constituents. 
Consistency requires that this evolution should account 
not only for the interparticle potentials but also for the 
effects of self-forces which encapsulate both conservative 
(backscattering of radiation from the system) as well as 
non-conservative (radiation recoil) effects. An effective 
field theory description of this procedure in the non- 
relativistic (or "post-Newtonian") limit was first formu- 
lated in [l|, and extended in to include higher order 
radiative corrections and in |5| to account for radiation 
reaction. 

Our focus in this paper is on the I = mass mode M 
and its time evolution due to gravitational radiation. Al- 
though this time dependence may be neglected to an ex- 
cellent approximation at low orders in the case of nearly- 
adiabatic, post-Newtonian systems, at higher orders the 
time dependence plays a role in the dynamics and must 
be consistently accounted for. 

We find that time dependence is necessary from the 
point of view of the renormalizability of the theory. Even 
if the mass parameter is made constant at some renormal- 
ization scale, renormalization group flow at the classical 
level induces time dependence. The fact that the mass 
of a radiating system changes with time implies that our 
mass is not the ADM mass. In fact, as we shall see, 
the boundary value of the mass, in terms of the renor- 
malization group scale, will be the Bondi mass. We also 
obtain a time evolution equation for the renormalized 
mass parameter. In order to do this, we impose conser- 
vation of a suitably defined energy-momentum (pseudo- 
censor operator for the composite system of gravitating 
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FIG. 1. Leading diagram topologies for the computation of 



sources, whose expectation value is computed using the 
closed time path formalism familiar from non-equilibrium 
quantum field theory and statistical mechanics |f|, here 
adapted to the classical system (see Q for a purely classi- 
cal formulation). Finally, as a byproduct of our analysis, 
we use the solution to the coupled renormalization group 
equations for M and to predict the pattern of loga- 
rithms that appear in the distribution of low-frequency 
I = 2 gravitons radiated from the bound system. 



RENORMALIZATION AND TIME 
DEPENDENCE 

In what follows, we work in the classical limit, tak- 
ing the moments in Eq. (JXJ) as c- number sources. Grav- 
itational dynamics is taken to follow from the Einstein- 



Hilbert action, Seh = ~ 16 ^ G J d 4 x^/gR. Our goal is 
to set up a causal time evolution equation for the point 
sources in Eq. ([I]) as they radiate gravitons. The cor- 
rect language for this sort of dissipative, causal prob- 
lem is the in-in formalism Q. We expand the metric 
SV = <?/it/ + h^y into a weak background g^ v = T]^ + h^ v 
and a fluctuating field and integrate out the fluctu- 
ations according to the in-in prescription. This yields an 
effective action 

e ir[h + ,h-,i,M] = J Dh + (x)Dh-(x)e iS ^ + ' M ^- iS ^~' M ' I \ 

(2) 

where we double the number of degrees of freedon and 
impose closed-time path boundary conditions on the fluc- 
tuating modes. In this paper we treat the multipole mo- 
ments as classical, so it is not necessary to double them 
for our purposes. The main object of interest is the real- 
time expectation value of the effective energy momentum 
tensor, which is 



(in\T^(x)\in) 



5T 



(3) 



h+=h-=0 



As a consequence of Ward identities, it obeys 
d l j,(T fJ, ' / (x)) = 0, which in fact fully determines the evo- 
lution of the moments. 

Diagrammatically, {T^J) consists of graphs with one 
external graviton (or equivalcntly h~). To first non- 
trivial order in the long wavelength approximation, the 



relevant diagram topologies are given in Fig. [T] Also 
at this order in the expansion, we may set the spatial 
momentum of the external graviton to zero in the calcu- 
lation of (T/iu). We also assume that the center of mass 
worldline is at rest at the origin x M (A) = (t, 0, 0, 0). Then 
the diagrams in Fig. Q] all give (T Ql ) of order V(S 3 (x), and 
from Fig. [2(a), 



(T 00 (x))=M(t)5 3 (x) 



(4) 



so that at this order the conservation law do(T 00 (x)) = 
simply yields M = 0. 

The diagrams in Fig. QJb), (c) have a more interest- 
ing structure. Although these diagrams do not contain 
logarithmic divergences in d — 4 spacctime dimensions, 
the computation of these terms does involve the one-loop 
integral 



id— i. 



(2tt) c 



r(3/2-d/2) „_ 3 
to 



(5) 



which is linearly ultraviolet divergent in d = 4. While 
such power divergences are set to zero within dimensional 
rcgularization in the MS scheme, other schemes which 
do not automatically set such divergences to zero would 
require the introduction of a time-dependent counterterm 
for the mass mode. Thus even if the time dependence of 
M is set to zero in the bare Lagrangian, renormalizability 
of the theory requires the mass to become dynamical. 

In the computation of the graphs in Fig. [Ijb), (c), it 
is sufficient to set h~ = 0, in which case one finds that 
the various propagators in the internal lines combine into 
retarded graviton propagators. Usingthe Feynman rules 
and the gauge-fixing choice in Ref. [1(, one finds the x- 
spacc result 



r(5)/ 



r(4)« r(l), 



df/W(f)/«(0). (6) 

-OO / 



G_ 

( c ) ~ 3Q " 



(7) 



Including the contribution of Eq. (|4]) this gives the time 
dependence of the I = mode 0, 



(T uu (x))=^(x) 



M 



G 



dfl^(f)l^(f) 



(8) 



The in-in expectation value is causal, depending only on 
the past history of the system. Energy-momentum con- 
servation at this order in the multipole expansion imposes 



1 For related results in a different formalism see 0, @ ■ 
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FIG. 2. Leading diagram topologies which yield logarithmi- 
cally UV divergent contributions to (T^). 



non-trivial time evolution for the mass mode, 



M(t) 



5 « 



m (t)I§\t) 



(9) 



As a check of our result, note that upon time-averaging 
both sides of this equation, we recover the textbook 
quadrupole radiation formula (AI) = — ^(-fy -fy ), re- 
flecting the fact that the time dependence of the mass 
accounts for energy loss into gravitational radiation. 



RENORMALIZATION GROUP FLOWS 

Ultraviolet divergences in four-dimensional gravity 
first appear in the order G 2 contributions to (T^ v ). The 
relevant diagrams are given in Fig. [2j We focus here on 
obtaining renormalization group equations, and therefore 
we only retain the parts of the diagrams that are singu- 
lar as d — > 4. Because the divergences of in-in correlators 
are the same as in the time-ordered scheme, it is sim- 
pler to employ time-ordered Feynman rules. Given the 
results of the previous section, the time dependence of 
M(t) would only yield contributions to the RG equation 
at order G 3 , and we may safely treat M as a constant 
in the calculation of the diagrams in Fig. @. Thus up 
to higher order corrections, we may replace M(t) by its 
time average M = J_ T M(t)dt/(2T) in the vertices. 

Working in momentum space with external graviton 
momentum U l = {ui, 0,0,0), we find that the divergent 
parts of the diagrams in Fig. [5] can be written as 



Fig. m- 



i-K 1 2 - [ dw 1 duj 2 

To G MI J 

x(27r)5(w-wi -LU 2 )P{LU 1 ,bJ 2 ) + ■■■ , (10) 



where we omit finite contributions. The function 
P(u>i, w 2 ) is a symmetric polynomial in its two arguments 



and X is the integral 
1 



L q d 



d-U 



1 



(27r) d " 1 (2 7 r) d - 1 q 2 p 2 (q + p) 



32tt 2 



d-4 



- In /i 2 + 



(11) 



with /1 an arbitrary subtraction scale. We find for each 
diagram 



P(wi,W 2 )(d) 
P(U1, W 2 )(e) 
P(W1,W 2 )(/) 



—32u> 2 ui 2 (w 2 + cjiu; 2 + wf) 1 
u)\ + u>l - 8(u>l u> 2 + OJ1LJ2) 
-49(cj 2 w 2 + u\u4) + 48w?w 2 : 



p(wi,w 2 ; 



(s) 



13 



83 



74 



+— (w 1 a;2 + wjw 2 ) + — Wi 



3, ,3 



(12) 
(13) 

(14) 
(15) 



P(wi,w a )(/,) = g(wi +w 2 ) 2 [-Mcj 2 oj 2 + 3(wJ + w£) 

+68(w 3 w 2 + wiw 3 )] , (16) 
giving the total result 

P(wi,w 2 ) = 32 [(wfw 2 + wicjg) + ^1^2 - (wiw 2 2 + WjWj)] ■ 

(17) 

The logarithmic divergence in the Feynman diagrams 
of Fig. [5] can be absorbed into the mass counterterm 
-£M(w) of Fig. [Ha), yielding a finite result for (T 00 ) 
at this order and a correponding renormalization group 
equation 



J ^^(^M 

ua)P(ui,U2), (18) 



G 2 M 
80 

x (27r)5(w — Wx 



or equivalcntly, the local x-space expression 



ft—M(t,ft) 
a/i 



2G M (iT&tM 
5 V" « 



or(4) r (2) 



r(3) r(3) 

(19) 



(*) 



In addition to mass renormalization, the theory de- 
fined by Eq. ([1]) also induces logarithmic renormalization 
of the quadrupole moment. This effect can be traced to 
singularities in graviton scattering off the gravitational 
field of the system, in particular from the (GM/r) 2 gen- 
eral relativistic correction to the potential. The result, 
first written in jij, is: 



(20) 



There is a simple interpretation of the mass RG flow. 
The mass as parameter as defined in Eq. ([1]) is the en- 
ergy of the conservative system not including energy in 
gravitational radiation, and it is not constant in time as 
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FIG. 3. An on shell graviton is emitted and scattered back 
( denoted by the small box) at a distance p. Observers at 
different distance would not agree on the value of the mass. 



we showed explicitly in the last section. At the order 
we are working it is clear that our mass definition is the 
Bondi mass, which gets no contribution from the radia- 
tion at infinity. That this mass is scale dependent can 
be roughly seen by noting that the logs in the mass arise 
as a consequence of radiation backscattering off of the 
static background. In Fig. [3] we illustrate such a phys- 
ical process. If the graviton scatters at a radius p then 
an observer at r < p would see radiation "at infinity", 
whereas an observer at R > p would include the energy of 
the backscattered graviton in the definition of the mass. 
Thus as observers reduce their distance to the center of 
the black hole system the mass parameter will diminish. 
In this sense the mass parameter is "asymptotically free" . 



APPLICATIONS 

For a gravitationally bound system of typical size 1Z 
we define a velocity parameter by v = TZui. For low fre- 
quency radiation, with r s u <»Cl where r s = 2GM 
is the gravitational radius, the RG running of Zy(u;,/i) 
is an 0[(r s u>) 2 } effect. On the other hand, the running 
of the mass mode is suppressed by v A relative to this ef- 
fect. Thus, to capture the leading logarithms of the form 
v 4 (r s ui) 2n In" v in the mass, we may solve the coupled 
RG equations by first inserting the solution to Eq. (|20| . 
assuming that M is /i-independent [H , 



Iij{(J,po), (21) 



where /?/ = —214/105 the beta function coefficient in 
Eq. (|20[) and Mq is the time-averaged mass, renormalized 
at the scale po. By time averaging Eq. (fT0|) . we find the 



running mass M(p) evolves according to the equation 



^lnM = -2G 2 </fif>„ 



(22) 



where is renormalized at the scale p. The 

solution is 



M(p) 
M 



exp 



//-(2)r(2)x _ /r (2) r (2)> 
Vij 1 ij /Mo Vij 1 ij In 



(23) 



where in terms of time averages at a high frequency scale 
Mo, 

/7-( 2 M 2 )\ _ >T ^ ir s> /r("+2) r ("+2)\ ln n M 

Vij 1 ij In — 2^ 2"n' [ ij ij >t "° u ' 



with r s = 2GMq. Taking p/po ~ TZlu = 
in the series are of the form [r s u>) In" v. 
leading logs in the running mass are 



(24) 

v, the terms 
The first few 



M(p) 

Mn 



1 



1 (/> ^ 



(3)r(3)v 



-ri lnw + 



107 (tf'L, 



(4) r(4)> 



420 



Ml 



11449 (Il-'l, 



(5) r (5)\ 



132300 



Ml 



-r\ In 3 v + 



In 2 v 



(25) 



and the term linear in the logarithm agrees with the re- 
sult of ll| obtained by different methods. 

This result can be used to determine logarithmic terms 
in the graviton emission amplitude from the system. In 
particular, the single (i = 2) emission amplitude to order 
r s u is given by (see e.g. Q) 



-4o(w) 



1 + 2-kGMuj, 



(26) 



where the leading order amplitude for polarized gravitons 
is A (oj) 



4m i 



-c*jlij(w, p), with tripi = 327rG. Using 



rcnormalization group invariance, we may evaluate the 
running quantities hj(y>) and M at a scale p ~ ui in or- 
der to minimize possible large logarithms in higher order 
contributions to |„4(w)| 2 . From the solution to the RG 



equation for the mass of Eq. ([23| , wc obtain 



A(u>) 



1 + irr s LU 



TT 

2 

(4) r(4h 



/r(3) r(3)v 

Vij 1 ij ;o 3 

-rzuj in v 



Ml 



W7tt (tt-'l, 



420 



Ml 



>° 5 , 2 , 

— r"u)m v + 



(27) 



The logarithmic corrections scale as irv A (r s uj) 2n+1 In" v, 
and we note that at this order, further logarithmic terms 
may arise from the subleading renormalization of the 
quadrupole moment. 

Unlike the first order contributions, the logarithmically 
divergent part of (T 00 ) at order G 2 is a conservative ef- 
fect. It therefore makes sense to include it in the defi- 
nition of a conserved energy E. For post-Newtonian in- 
spiralling binary systems with quasi-circular orbits, the 



5 



conserved energy is commonly expressed as a function of 
the orbital frequency f2. Logarithms in E(Q) arise first at 



fourth post-Newtonian order [11[ . There is a direct con- 
tributions from the energy we computed here as well as an 
indirect contribution due to the the use of the equations 
of motion in deriving the relation between orbital radius 
and frequency. The equations of motion may be inferred 
from energy conservation, by essentially reversing the ar- 
guments in [ll|. We find the leading contribution linear 
in logarithms to be 



E(Cl) 



H 448 
2 ~15~ 



vx In x 



(28) 



where /i is the reduced mass, v = h/Mq and x = 
(GM n) 2 / 3 . This is in agreement with the 4PN loga- 
rithm computed in [ill ]. We may also use our result for 
the running mass in Eq. f|23[) or ([25)) to extract higher 
powers of the leading logarithms in E(Q) by resumming 
the leading logarithms to all orders [12j |. which includes 
the In™ x terms at (4 + in) post-Newtonian order includ- 
ing their numerical coefficient. 



CONCLUSIONS 

In this paper, we have set up a formalism for determin- 
ing the real-time evolution of the dynamical moments 
that describe a system of gravitationally bound black 
holes. In order to account for the dissipation of energy 
due to radiation in a way that respects causality, the in-in 
formulation of quantum field theory must be employed. 
We find that the time evolution of the i = mass mode is 
inextricably linked to the renormalization of the effective 
theory. At second order in the gravitational coupling G 
one finds logarithmic ultraviolet divergences even in the 
classical theory, which induce non-trivial RG flows. The 
RG equation for the mass mode, together with previously 
obtained results for the running of the I = 2 moment, re- 
sums the leading logarithms in the conservative energy 
of the form v i (r s uj) 2n In™ v. 

The results of this paper can be extended in several 
of directions. It is clear that the methods introduced 
here can also be used to set up evolutions equations for 
the higher moments (the center of mass momentum and 



angular momentum as well as higher moments) that take 
into account radiative losses. To do so would require 
computing the similar diagrams to those in Fig. [TJ but 
with non-zero spatial momentum. 

Another direction of research would be to consider the 
evolution of quantum black holes within this formalism, 
and the consequences of the RG flows for (e.g. Hawking) 
radiative processes. In that case, we expect that at lead- 
ing order in the multipole expansion, the corresponding 
evolution equation for the mass operator will take the 
form of a relation between its expectation value and cer- 
tain two-point correlators of the quadrupole moments, 
some of which have already been obtained in [2[ in the 
low-frequency limit. 
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(AR, IZR). 
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